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Systematic investigation of the accuracy of the description of the energies of deformed one-quasiparticle states has been performed 
in covariant density functional theory in actinide and rare-earth mass regions. The sources of the discrepancies between theory 
and experiment are analyzed. Although some improvements in the description of ground state configurations and one-quasiparticle 
spectra can be achieved by better parametrization of the relativistic mean field Lagrangian, the analysis suggests that spectroscopic 
fxj quality of their description can be achieved only in theoretical framework which takes into account particle-vibration coupling. 
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1. Introduction 

Further progress in understanding low energy nuclear phe- 
nomena in stable and exotic nuclei is strongly connected with 
the development of nuclear density functional theory (DFT) 
in its non-relativistic and relativistic (covariant) incarnations. 
These theories provide rather successful description of differ- 
ent aspects, such as deformations, masses, collective excitations 
etc, of nuclear many-body problem, see Refs. AH 

and refer- 
ences therein. In addition, they are aimed on global (i.e. across 
full nuclear chart) description of the nuclear properties. 

However, absolute majority of the applications of nuclear 
DFT has been focused on collective properties of nuclei. This is 
due to the fact that time-odd mean fields are needed for the de- 
scription of one-(many)-particle configurations which are char- 
acterized by broken time-reversal symmetry in intrinsic frame 
(see Ref. |2] and references therein). As a result, the descrip- 
tion of such configurations is more complicated as compared 
with the one of the ground states of even-even nuclei. 

There are only few features of deformed nuclear systems, 
strongly dependent on single-particle degrees of freedom, 
which have been addressed in the DFT studies on the mean- 
field level and compared with experiment. When discussing the 
impact of single-particle degrees of freedom on nuclear proper- 
ties, it is important to separate different aspects of their physics 
which only weakly depend on each other. These are (i) de- 
formation polarization effects induced by particle or hole, (ii) 
alignment properties of single-particle orbital in rotating poten- 
tial and relevant polarization effects in time-odd mean fields, 
and (iii) the single-particle energies. 

Let us consider each of those aspects. The addition of the 
particle to (or the creation of the hole in) even-even core in- 
duces deformation polarization effects. The investigation of su- 
perdeformed rotational bands in the A ~ 140 - 150 mass re- 
gion, which are one of the best examples of undisturbed single- 
particle motion Si , has revealed that deformation polariza- 



tion effects induced by particle or hole are well described in 
nonrelativistic S \^ and covariant S 01 DFT's; the average 
deviation from available experimental differential quadrupole 
moments is around 20%. Similar (but somewhat less accurate 
because of the role of pairing) results have been obtained also in 
the A ~ 130 mass region of high- and superdeformation Q]. 
Macroscopic+microscopic (MM) method based on the Nilsson 
potential describes deformation polarization effects reasonably 
well but suffers from the fact that these effects are not uniquely 
defined fli]. 

Alignment properties of single-particle orbital in rotating po- 
tential can be accessed via effective (relative) alignments 11011 
of two compared rotational bands. This quantity sensitively de- 
pends on both the alignment properties of single-particle or- 
bital by which two bands differ and polarization effects (mostly 
in time-odd mean fields) induced by the particle in this or- 
bital II ill . Effective alignments are on average better repro- 
duced in the covariant DFT (CDFT) calculations than in the 
cranked Nilsson-Strutinsky version of the MM approach based 
on phenomenological Nilsson potential, see comparisons pre- 
sented in Refs. jl llll2lll3ll . Reasonable description of effective 
alignments can be obtained also in Skyrme DFT S Q3], but 
it is somewhat plagued by the uncertainties in the definition of 
the coupling constants for time-odd mean fields lfl4l [l5ll . This 
problem does not exists in the CDFT since time-odd mean fields 
are defined via Lorentz covariance H. 

The discussion above clearly shows that some aspects of the 
single-particle motion are described on average better and in a 
more consistent way in self-consistent DFT than in phenomeno- 
logical MM method. This is despite the fact that no single- 
particle information (apart of some spin-orbit splittings in the 
case of Skyrme DFT) has been used in the fit of the DFT param- 
eters, while the parameters of the Nilsson or Woods-Saxon po- 
tentials are fitted to experimental sing le-particle energies (see, 
for example, Refs. IU6lll7lll8lll9l l201). However, starting from 
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Figure 1: The evolution of one-quasineutron energies as a function of neutron number for the Pu isotopes. Hole states are plotted below the ground state (zero 
energy), and particle states are plotted above. Experimental data (one-quasineutron band-head energies) are taken from Ref. 13211 . The squared amplitudes of the 
dominant component of the wave function change gradually with neutron number. Thus, their variations are presented below in the neutron number range where 
these states are plotted in left panel in the following format: [state label: A-B], where A is the squared amplitude of the dominant component of wave function at 
lower value of N and B at upper value of N. These variations are 1 /2[501] : 0.67 - 0.46, 3/2[501] : 0.82 - 0.87, 5/2[752] : 0.57 - 0.56, 5/2[633] : 0.77 - 0.78, 
3/2[631] : 0.55-0.57, 7/2[743] : 0.72-0.74, l/2[632] : 0.55 - 0.57, 7/2[743] : 0.72-0.74, 1/2[631] : 0.58 - 0.51, 5/2[622] : 0.65 - 0.69, 7/2[624] : 0.86-0.89, 
7/2[613] : 0.78 - 0.83, 3/2[622] : 0.66 - 0.52, 9/2[615] : 0.96 - 0.96, 9/2[734] : 0.82 - 0.85. The l/2[7] state is strongly mixed. However, the cumulative squared 
amplitude of the components of the wave function with N = 7 in the structure of this state exceeds 90%. Thus, we label it only by principal quantum number N and 

n. 



earlier studies of the single-particle spectra of few nuclei in ac- 
tinides within covariant [21] and Skyrme DFT B22I1 . and follow- 
ing by later global survey of the ground state configurations in 
odd-mass nuclei (Ref. B23IH and the investigations of the spectra 
of odd-proton Ho nuclei (Fig. 6 in Ref J150) in triaxial Skyrme 
DFT and the spectra of selected Rb lET Y and Nb JUj nu- 
clei in axial Gogny DFT i25ll . it became clear that the single- 
particle spectra are poorly described with modem DFT. It is 
obvious that these theories do not possess spectroscopic qual- 
ity description of the single-particle spectra which is achievable 
in the MM method as a consequence of the fit of the potential 
parameters to experimental single-particle energies. 

Despite these results, no systematic analysis of the accuracy 
of the description of the single-particle spectra has been per- 
formed in DFT. The current manuscript aims on such analysis 
within the CD FT using extensive set of experimental data on the 
energies of one-quasiparticle states in deformed nuclei. There 
are two main goals behind of this study. First, typical uncertain- 
ties in the description of single-particle energies will be defined. 
Second, the sources of discrepancies and possible ways of the 
improvement of the description of the single-particle spectra 
will be discussed. 



2. The details of calculations 

Cranked relativistic Hartree-Bogoliubov (CRHB) approach 
is used in the current manuscript in a similar way as it was done 
in an earlier study of the spectra of few odd-mass actinide nu- 
clei in Ref. 12111 . Time-odd mean fields are taken into account, 
and the blocking due to odd particle (see Sec. VI. A in Ref. Q21 ] 
for details) is performed in a fully self-consistent way. The D 1 S 
Gogny force is used in the pairing channel. Nuclear configura- 
tions of deformed odd nuclei (further one-quasiparticle [1-qp] 
configurations) are labelled by means of the asymptotic quan- 
tum number Q[Nn z A] (Nilsson quantum number) of the domi- 
nant component of the wave function of blocked single-particle 
orbital. In each nucleus under study, the binding energies are 
calculated for a number of the 1-qp configurations based on the 
orbitals active in the vicinity of the Fermi level, and then the 
1-qp spectra are built as shown in Fig. Q] 

The CRHB equations are solved in the basis of an anisotropic 
three-dimensional harmonic oscillator in Cartesian coordinates. 
The same basis deformation /3o = 0.3, y — 0° and oscillator 
frequency huto = 41A -1 ^ 3 MeV have been used for all nuclei. 
All fermionic and bosonic states belonging to the shells up to 
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Figure 2: The distribution of the deviations of the calculated energies 
E„p(theor) of one-quasiparticle states from experimental ones E qp (exp). The 
vertical axis shows the percentage of the states which deviate from experiment 
by the energy deviation range (the width of bar) specified on horizontal axis. 



Nf = 14 and Nb = 18 are taken into account in the diago- 
nalization of the Dirac equation and the matrix inversion of the 
Klein-Gordon equations, respectively. The comparison with the 
results of calculations in a larger fermionic basis (Nf = 16) 
shows that the energies of quasiparticle states are described 
with an accuracy of approximately 50 keV which is sufficient 
for a statistical analysis. 

The calculations have been performed with the NL1 ll26ll and 
NL3* 12711 parametrizations of the relativistic mean field (RMF) 
Lagrangian. The NL1 parametrization has been fitted to the nu- 
clei in the valley of beta-stability. On the contrary, the fit of 
the NL3* parametrization includes neutron-rich nuclei so it is 
partially tailored towards the description of such nuclei. This 
recently fitted parametrization has been successfully applied to 
the description of binding energies 112711 . ground state proper- 
ties of deformed nuclei [28], spectra of odd spherical nuclei 
within the relativistic particle-vibration model 113011 ^ rotating nu- 
clei [27], giant resonances 12711 . and breathing mode [31]. Note 
that only bulk properties of nuclei such as binding energies, 
radii etc. and nuclear matter properties have been used in the fit 
of these two parametrizations. 



Table 1: The summary of calculations. The number of calculated 1-qp configu- 
rations and the number of experimental 1-qp states compared with calculations 
are displayed in second and third columns, respectively. Fourth column gives 
the percentage of ground states, the structure of which is correctly reproduced 
in the calculations. 
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Figure 3: The same as Fig. |2] but for the case when the energy scale of theoret- 
ical spectra is corrected for low Lorentz effective mass. 



nuclei with octupole deformation have been excluded from con- 
sideration. In addition, the calculations have been performed 
in rare-earth region. However, the experimental data on 1-qp 
states in rare-earth region are significantly larger than in ac- 
tinide region 0211 . As a consequence, the calculations were 
performed only for odd-proton Z = 61 (Pm), 63 (Eu), 65 (Tb), 
67 (Ho) and 69 (Tm) isotope chains and only with the NL1 
parametrization of the RMF Lagrangian; in total, for 31 odd- 
proton nuclei. All these nuclei are deformed with quadrupole 
deformation fi^ > 0.2 and labelling of their single-particle states 
by means of the Nilsson quantum numbers is commonly ac- 
cepted [133, 34, 35, 36, 37, 38]. Table 1 provides the summary 
of these calculations. The data on adopted experimental one- 
quasiparticle levels are taken from Ref. [32]. 



1-qp spectra were calculated for 44 nuclei, namely, 21 odd-Af 
and 23 odd-Z nuclei in actinide region with Z = 89 - 100. The 



3. Statistical analysis 

Fig. Q] shows the comparison of calculated and experimental 
one-quasineutron spectra in Pu isotopes. A number of features 
are clearly seen. First, for a given 1-qp state the discrepancy be- 
tween theory and experiment depends on neutron number. Sec- 
ond, for a given 1-qp state the slope of the energy curve as a 
function of neutron number is more pronounced in the calcula- 
tions than in experiment. These two features are interconnected 
and they emerge from the fact that theoretical energy scale is 
more stretched out than experimental one due to low effective 
mass (see discussion below). The change of the Fermi energy 
with neutron number leads to the changes of the energy differ- 
ences between ground and excited states and these differences 
are affected by effective mass in the calculations. Third, the 
relative energies of different experimental 1-qp states are not 
always reproduced in calculations. This feature originates from 
the fact that the energies of spherical subshells, from which 
deformed states emerge, can deviate from experiment ll2l1l ). 
These three features are seen in all isotone and isotope chains. 

Note however that for a given state the 1-qp energy changes 
as a function of particle number are appreciably smaller when 
1-qp energies of proton (neutron) subsystem are shown as a 
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function of neutron (proton) number (see, for example, Fig. 6 
in Ref. 11511 ) because the changes in proton (neutron) Fermi en- 
ergy and deformation are relatively small. 

Statistical analysis of the discrepancies between calculated 
and experimental energies of one-quasiparticle states is pre- 
sented in Fig. 12 One can see that in the actinide region only ap- 
proximately 33% of one-quasiparticle states are described with 
an accuracy better than 200 keV, and approximately 22% with 
an accuracy between 200 and 400 keV. The percentage of states 
for a given range of deviations goes gradually down with an in- 
crease of the deviation between experiment and calculations. 
However, for some states the deviation of calculated energy 
from experimental one exceeds 1 MeV and can be close to 1 .4 
MeV. Note that earlier analysis of the single-particle spectra of 
few actinide nuclei in Ref. J2~lll shows similar features. Fig. [2] 
also shows that with the NL1 parametrization the 1-qp energies 
in odd-proton rare-earth nuclei are somewhat better described 
as compared with actinide region. For example, the energies of 
44% of the states are described with an accuracy better than 200 
ke V. Otherwise, the distrubion histograms for the deviations are 
similar in both regions and for both parametrizations. 

Fourth column of Table[T]shows the percentage of the ground 
states, the structure of which is correctly reproduced in the cal- 
culations. These values are comparable or somewhat better 
than the ones obtained in systematic Hartree-Fock+BCS calcu- 
lations of deformed nuclei with SHI, SkM* and SLy5 Skyrme 
forces and FRDM calculations employing phenomenological 
folded- Yukawa potential [23], which show the agreement with 
experiment at about 40% level. 

It is interesting that the overall accuracy of the description 
of the energies of deformed one-quasiparticle states is slightly 
better in old NL1 parametrization, which was fitted 25 years 
ago mostly to the nuclei at the y3— stability line, than in the re- 
cent NL3* parametrization. This suggests that the inclusion of 
extra information on neutron rich nuclei into the fit of the NL3* 
parametrization may lead to some degradation of the descrip- 
tion of single-particle states along the valley of beta-stability. 



4. The sources of the discrepancies between theory and ex- 
periment 

4.1. Effective mass 

Low effective mass of CD FT is one of the sources of the dis- 
crepancies between theory and experiment. This is because the 
average level density of the single-particle states on the mean 
field level is related to the effective mass m* L (kf)/m (Lorentz 
mass in the notation of Ref. [39] for the case of CDFT) of the 
nucleons at the Fermi surface which depends on momentum as 



m* (k F ) 



(m*(0) 



(m*(0) 



+ 0.08 



(1) 



where m*(Q)/m is the value at k — which is called Dirac ef- 
fective mass fl39ll and 0.08 is obtained for {kf/m) 2 when typical 
value kf « 1 .35/fm is used. The values of Lorentz effective 
mass for employed parametrizations are given in Fig. [3] 

Low effective mass leads to a stretching of theoretical single- 
particle energy scale as compared with experiment (see Refs. 
B29L BOIl for comparisons of calculated and experimental spec- 
tra in spherical nuclei). The role of the energy stretching due 
to low effective mass can be illustrated on the example of the 
7rl/2[420],7r3/2[411] and;r5/2[402] states in rare-earth region. 
These states emerge from the nds/2 spherical subshell. For the 
majority of the nuclei understudy, the 7r3/2[41 1] state is located 
close to the proton Fermi level. As a consequence, it is observed 
in 25 nuclei, and its energy in those nuclei is described in the 
calculations with an average accuracy of 250 keV. On the con- 
trary, the 7rl/2[420] (7t5/2[402]) state is located significantly 
below (above) the Fermi level both in experiment and calcu- 
lations. However, stretching of theoretical energy scale due to 
low effective mass results in systematic deviations of calculated 
energies of these states from experiment by more than 1 MeV. 
This leads to a peak in distribution histogram of the deviations 
of calculated energies from experimental ones at the deviation 
energy of around 1.1 MeV (see Fig. |2j. Another contributor to 
this peak is the 7r9/2[514] state. 

This stretching is also clearly visible when single-particle 
Nilsson diagrams obtained in the CDFT and phenomenological 
Nilsson or Woods-Saxon potentials are compared (see, for ex- 
ample, Fig. 1 in Ref. 11211 '). These potentials are characterized 
by an effective mass m*(kf)/m sa 1.0 which gives calculated 
level density close to experiment. As illustrated in spherical 
nuclei, the calculated level density and single-particle spectra 
are compressed and come closer to experimental ones when 
CDFT is suplemented by particle-vibrational coupling (PVC) 
11291 13011 (similar effect exists also in the PVC models based on 
non-relativistic DFT HiHlil). 

Similar compression of calculated spectra is expected also 
in deformed nuclei. However, no PVC model based on DFT 
framework has been developed so far for such nuclei. System- 
atic comparison of single-particle Nilsson diagrams obtained in 
the CDFT and phenomenological Nilsson potential suggest that 
on average expected compression of the single-particle spectra 
can be achieved by multiplying the energies of the 1-qp con- 
figurations by Lorentz effective mass. Such energy rescaling 
follows also from the results of phenomenological scheme of 
Ref. [43] based on a linear ansatz for the energy dependence 
of the scalar and vector components of the nucleon self-energy 
for the states close to the Fermi surface which simulates the dy- 
namical effects that arise from the coupling of single-nucleon 
motion to collective surface vibrations. 

The impact of such energy rescaling on the distribution of 
the deviations between theory and experiment is shown in Fig. 
[3] One can see that more than 75% of states are described with 
an accuracy better than 400 keV; this is a typical accuracy of the 
description of the energies of the deformed 1-qp states within 
phenomenological potentials [36, 44]. Although this energy 
rescaling is somewhat schematic, it clearly indicates that PVC, 



leading to an increase of effective mass 



irly it 
[29], 



will also improve 



the description of experimental spectra as compared with mean 
field results; this has already been seen in spherical nuclei 113011 . 



5. Consequences for spectroscopic quality energy density 
functional 



4.2. Relative energies of different states 



Incorrectly calculated relative energies of different states (see 
Fig. [TJ represent another source of the discrepancies between 
theory and experiment (see discussion of the Pu isotope chain 
above). It originates from the fact that relative energies of dif- 
ferent spherical subshells, from which deformed states emerge, 
are not properly reproduced in model calculations lEfl . This 
source shows up also in a statistical analysis. For example, the 
energy of the hole vl /2[501] state in actinide region is system- 
atically higher than experimental one by around 1 MeV in both 
employed parametrizations. This is only deformed state origi- 
nating from the vpx/2 spherical subshell, so the current analysis 
suggests that in the calculations it is placed too deep (by ap- 
proximately 1 MeV) with respect of other spherical subshells. 
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Figure 4: Histogram of differences between the moment of the i-th 1-qp con- 
figuration and the average moment of all calculated 1-qp configurations in the 
(N,Z) nucleus. The distributions for mass quadrapole (panel (a)) and mass 
hexadecapole (panel (b)) moments are shown. 



The need for quantitative understanding and highly accurate 
description of nuclear structure phenomena is a driving force 
behind the efforts to develop spectroscopic quality energy den- 
sity functionals (EDF). Such efforts are especially visible in the 
field of non-relativistic Skyrme DFT IU5ll46ll . However, the ba- 
sic question whether this is possible on the mean-field level is 
still under debate [45]. 

It is well known from Skyrme DFT that EDF, forced to 
describe accurately single-particle spectra on the mean-field 
level, are characterized by effective mass close to one 14711 . In 
the CDFT, one cannot improve the description of the single- 
particle spectra by increasing effective mass since all CDFT 
parametrizations on the Hartree level have Lorentz effective 
mass m* L {kF)lm close to 0.65 J2l. The current analysis for 
deformed nuclei strongly suggests that further progress in im- 
proving spectroscopic quality of covariant EDF will be quite 
limited on the mean-field level. We expect that even for other 
modern CDFT parametrizations, not employed in the current 
manuscript, the distributions of the deviations of calculated en- 
ergies from experimental ones of the 1-qp states will be compa- 
rable or only slighthtly better than the one seen in Fig. |2] This 
is due to their low Lorentz effective masses and unavoidable 
errors in relative placement of specific single-particle orbitals. 
However, for some parametrizations such as NLSH and NL- 
RA1 obtained distributions can be even worse than that of Fig. 
|2j this follows from the analysis of Ref. J2lll . 

As a consequence, the only way to substantially improve the 
description of the single-particle properties in the framework 
based on CDFT is to take into account PVC It was already 
illustrated in spherical nuclei that this leads to a significant im- 
provement in the description of the experimental energies of 
the dominant single-particle states [29, 30]. In addition, it takes 
care of well-known fragmentation of the single-particle strength 
of the levels; this feature is completely ignored on the mean- 
field level. Based on experience in spherical nuclei, one can 
expect that the inclusion of PVC in deformed nuclei will bring 
Lorentz effective mass close to one, thus leading to a level den- 
sity which is similar to experimental one. This will definitely 
improve the description of the 1-qp spectra. The corrections to 
the energies of single-particle states due to PVC are strongly 
state-dependent in spherical nuclei 1I29I l30ll . On the contrary, 
in deformed odd nuclei the corrections to the energies of 1-qp 
states due to PVC are expected to be less state-dependent since 
surface vibrations are more fragmented in deformed nuclei than 
in spherical ones [ 38 HI- However, they can still affect the rela- 
tive energies of different deformed 1-qp states. 

6. Deformation effects 

Each blocked single-particle orbital induces deformation po- 
larization effects. Fig. [3] shows the distribution of the mass 



This provides extra justification for the energy scaling procedure employed 
in Sect. RTTl 
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quadrupole and hexadecapole moments of the 1-qp configura- 
tions relative to the respective average moments of the full set 
of 1-qp configurations in a given nucleus. This figure clearly 
shows that actinide nuclei are reasonably rigid since polariza- 
tion effects (especially for quadrupole moments) induced by 
odd particle do not lead to substantial deviations of equilibrium 
moments from average values. On the contrary, rare-earth nu- 
clei are considerably softer than actinides since their histograms 
are significantly wider. The results for actinides show also that 
polarization effects do almost not depend on the parametriza- 
tion of the RMF Lagrangian; small differences between his- 
tograms obtained in the NL1 and NL3* parametrizations are 
most likely due to slightly different sets of the 1-qp configura- 
tions obtained in the calculations. These changes in the mo- 
ments/deformations and relevant changes in the binding ener- 
gies of the 1-qp configurations induced by odd particle have 
to be taken into account when experimental data are compared 
with the results of the calculations. However, they are com- 
pletely ignored in the models most frequently used for the anal- 
ysis of the spectra of deformed odd nuclei such as MM model 
1 44], particle +rotor [37] and quasiparticle-phonon ll38ll models; 
these models assume the same deformation for all 1-qp states 
in the nucleus under study. 

The calculations suggest that the deformation driving effects 
induced by odd proton or neutron are sufficient to polarize the 
nucleus towards positive or negative y-deformation in some 
one-quasiparticle configurations. On average, only 6.8% of cal- 
culated states have y-deformation in the range 1° < \y\ < 17°; 
the -/-deformation of other states is either exactly zero (in the 
majority of the cases) or below 1°. It turns out that the appear- 
ance of sizable y-deformation correlates with blocking of few 
specific single-particle states. For example, in the actinide re- 
gion 43% of triaxial configurations has the dominant v[622]3/2 
structure, and 38% has the v[501]l/2 structure. Time-odd mean 
fields have almost no effect on the equilibrium deformations 
similar to previous CDFT studies. On the contrary, the results 
of Skyrme EDF studies of Ref. Il5ll show that the inclusion of 
time-odd mean fields favors the axial deformation of calculated 
configurations. 

7. Conclusions 

In conclusion, for the first time a systematic analysis of the 
accuracy of the description of the energies of one-quasiparticle 
configurations in deformed odd nuclei has been performed in 
the DFT framework. It provides theoretical estimates on the er- 
rors in calculated energies of one-quasiparticle configurations. 
Two sources of inaccuracies, namely, low effective mass lead- 
ing to a stretching of the energy scale and incorrect relative po- 
sitions of some single-particle states exist in model calculations. 
While the solution of the latter problem can be attempted in the 
DFT framework, we do not believe that it will lead to signif- 
icant improvements in spectroscopic quality of energy density 
functionals. The comprehensive solution requires taking into 
account particle-vibration coupling which will (i) compress the 
calculated one-quasiparticle spectra bringing them closer to ex- 
periment and (ii) correct the relative energies of different single- 



particle states. In our opinion, only such approach combined 
with respective re-parametrization of the RMF Lagrangian can 
lead to spectroscopic quality energy density functionals. 

In addition, our results show that one should be extremely 
careful in the interpretation (predictions) of (for) the data which 
involves absolute or relative energies of different one-/two- 
/many-quasiparticle states since their energy description is as- 
sociated with non-negligible theoretical errors. Extra care has 
to be taken also in the case when physical phenomenon sen- 
sitively depends on the single-particle energies. For example, 
the existence of very shallow left and right chiral minima with 
depth of around 200 keV separated by a tiny barrier ( » 60 keV) 
is responsible for chiral rotation [14811 . Theoretical errors in the 
description of single-particle energies can either create false 
minima (barrier) or kill real ones. 
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